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STRESSES IN A TWO-BAY NONCIRCULAR CYLINDER UNDER TRANSVERSE LOADS' 


By George E. Griffith 


SUMMARY 

A method, taking into account the effects oj flexibility and 
based on a general eighth-order differential equation, is ■pre- 
sented for finding the stresses in a two-bay, noncircular cylinder 
the cross section of which can be composed of circular arcs. 
Numerical examples are given for two cases of ring flexibility 
for a cylinder of doubly symmetrical ( essentially elliptic) cross 
section, subjected to concentrated radial, moment, and tan- 
gential loads. The results parallel those already obtained for 
shells with circular rings. 

INTRODUCTION 

In airplane fuselages with flexible rings subjected to con- 
centrated loads, the stresses in the neighborhood of the load 
differ markedly from those given by the simple engineering 
formulas, and more refined methods, which take into account 
the interaction of rings and shell, are needed to predict the 
stresses accurately. The first paper on this subject, pub- 
lished in May 1944, was that of Wignot, Combs, and Ensrud 
(reference 1), who treated the circular cylinder subjected to 
concentrated loads but neglected the effect of the extensional 
deformations of the shell. Hoff (reference 2) gave a more 
complete analysis, including the effects of many rings, for 
the case of symmetric transverse loads. The results were 
corroborated experimentally by Kuhn, Duberg, and Griffith 
(reference 3), who also extended the theory to include con- 
centrated moment and tangential loads. Later, Duberg and 
Kempner (references 4 and 5) reduced the labor of computa- 
tion by giving the results in the form of charts and showed 
that for practical purposes it was usually sufficient to con- 
sider only a region within 2 bay lengths of the load. Further 
investigations considered additional effects, heretofore neg- 
lected, such as the shearing and axial deformation of the rings 
(references 6 to 8), bending rigidity of stringers (references 
8 and 9), shear carried by the stringers (reference 7), and 
eccentricity of ring and sheet (references 6, 7, 8, and 10). 

All investigations referred to dealt exclusively with rein- 
forced monocoque cylinders of circular cross section. The 
present report gives an analysis for a two-bay noncircular 
cylinder, enclosed between ring bulkheads rigid in then- 
planes, with the middle, flexible ring subjected to concen- 
trated and distributed loads. The fundamental assumptions 
used in the analysis are the same as those previously used fox- 
circular cylinders. 


Many noncircular fuselages can be closely approximated 
by using circular sections of different radii and joining the 
sections at points of tangency. The rings discussed herein 
are of this form. Associated with each i-ing section is a two- 
bay panel (fig. 1), any number of similar panels composing 
the complete two-bay cylinder. The solution for the stresses 
in such a structure is based on the development of a general 
eighth-order differential equation, wi-itten in tex-ms of the 
moment at the skin center line. A separate differential 
equation of the same general fox-m applies to each curved 
panel in the structure. Application of the correct boundary 
conditions results in sets of simultaneous equations which 
yield the unknown constants in the moment expression. 

Although a two-bay cylinder does not conform to the usual 
fuselage structure, it is believed, on the basis of comparisons 
with some of the work previously mentioned, that the results 
obtained ai’e indicative of those found in more complicated 
structures. In accordance with the findings of prior investi- 
gations, shear and axial deformations of the loaded x-ing are 
neglected, but eccentricity of ring and sheet is included. 

The numerical examples deal with a doubly symmetrical, 
two-bay cylinder of nearly elliptic cross section (fig. 2), sub- 
jected to concentrated loads. 

SYMBOLS 

„ l'R 6 

A =m 

Et'R 2 
1 GtU 

C, coefficient of shear flow in sheet 

C H coefficient of ring axial force 

Cm coefficient of ring bending moment 

C r coefficient of ring shear force 

E Young’s modulus in sheet and ring, pounds per 

square inch 



Figure 1. — Two-bay panel showing positive directions of displacements and sheet stresses. 


1 Supersedes NACA TN 2512, “Stresses in a Two-Bay Noncircular Cylinder Under Transverse Loads" by George E. Griffith, 1951. 
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Fhujre 2. — Two-bay cylinder, composed of four circular panels, with doubly symmetrical 
noncircular rings, used for numerical examples. 

G shear modulus in sheet, pounds per square inch 

H axial force acting on sheet, pounds 

1 moment of inertia of ring section, inches 4 

L length of bay, inches 

M bending moment at sheet center line, inch-pounds 
M c applied concentrated moment, inch-pounds 
M s bending moment about neutral axis of ring, inch- 
pounds 

P applied concentrated radial load, pounds 

R radius to sheet center line of circular panel, inches 

R i, 7?o radii of circular panels composing cylinder, inches 
T applied concentrated tangential load, pounds 

V shear force, pounds 

a±bi complex roots of auxiliary algebraic equation (used 
when real roots also occur) 

Ci±(hi, complex roots of auxiliary algebraic equation (used 
c 2 ±c/ 2 i when no real roots occur) 

e eccentricity of ring and sheet (distance between 

sheet center line and neutral axis of ring), inches 
j distributed radial load acting on sheet, pounds per 

inch 

h distributed axial load acting on sheet, pounds per 

inch 

;=v-t 

k lt k 2 real roots of auxiliary algebraic equation 
m distributed moment acting on sheet, inch-pounds 

per inch 

p normal force in stringers, pounds per inch 

2 shear flow in sheet, pounds per inch 

s distance in circumferential direction, inches 

t thickness of sheet, inches 

t' thickness of all material carrying bending stresses 

in panel if uniformly spread around circum- 
ference, inches 

u displacement of sheet in axial direction, inches 


v displacement of sheet in circumferential direction, 

inches 

v Q displacement of sheet in circumferential direction 

at ring, inches 

w displacement of sheet or ring in radial direction, 

inches 

x distance in axial direction measured from loaded 

ring, inches 

y shear strain in sheet 

e s normal strain in sheet in circumferential direction 

e z normal strain in sheet in ^-direction 

9 angular distance, degrees or radians 

a stringer normal stress, pounds per square inch 

cr s normal stress in sheet in circumferential direction, 

pounds per square inch 

r shear stress in sheet, pounds per square inch 


GENERAL DIFFERENTIAL EQUATION 

The basic element of the present analysis is a two-bay 
panel, as shown in figure 1, composed of sheet, longitudinal 
stiffeners or stringers, and transverse stiffening ring sections 
of constant radius. Distributed or concentrated loads are 
applied in the plane of the middle or flexible ring. By 
joining several panels at points of tangency, many two-bay 
cylinders of various cross-sectional shapes can be achieved. 

In actual practice the sheet covering is outside the rings, 
and this eccentricity of sheet and ring is henceforth taken 
into account. For convenience the moment considered is 
the moment which exists at the sheet center line rather than 
at the neutral axis of the ring. From this moment the 
bending moment in the ring is easily found. 

The underlying assumptions used in the analysis are as 
follows: 

(1) The shear stress, carried by the sheet alone, may vary 
in the circumferential direction but remains constant in the 
axial direction. 

(2) The material in the cross section of the panel (sheet 
and stringers) capable of carrying normal stresses due to 
bending of the panel is assumed spread around the circum- 
ference in a fictitious sheet of thickness l'. 

(3) The loaded ring has no torsional stiffness or bending 
stiffness out of its plane. 

(4) The end ring supports are restrained from deforming 
in their planes but are free to warp out of their planes. 

Under these assumptions, for any panel with constant 
geometrical properties, a general differential equation is 
developed for the moment at the sheet center line. All 
forces, stresses, and displacements in the panel may be 
obtained from that moment. The general differential equa- 
tion (see the appendix for development) is 


tPM 
de 8 ' 




3 B- 


-(2 — 3 B) 
2eA/ 


R 


d 6 M ( _ 

c w 1 v 

-6B+ 

6e 2 A.\ 
R 2 ) 

d 4 M 

de i 

_ e , e 2 \ 

ld 2 M 

+ 6^1| 

f e 

2 ~~R + R 2 ) 

J de 2 

1 — — 
V R 


I M-- 


■- F(e ) 
(i) 
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where 


SOLUTION OF GENERAL DIFFERENTIAL EQUATION 


A»)-« ! [^+a-3S)^-3S§]+ 

8 !S+( ! - !S i»+( 1 - m +t)» 


When the ratio e/R, a measure of the eccentricity of the 


ring and sheet, is very small 


(WD 


it can be neglected; 


equation (1) then contains only the two nondimensional 
parameters A and B, which are related to the geometry and 

t'R e Et'R 2 

physical properties of the panel: A=j^ 3 and B= CtL 2 ' 


The ratio A/B is an index of the over-all flexibility of the 
structure. When this ratio is very small the ring section is 
rather rigid. A large value of A/B (500, for example) indi- 
cates that the ring section is somewhat flexible, so that radical 
departures of the stress distributions from the elementary 
values result. In actual practice the parameter B varies 
over a small range (between about 10 and 80), whereas A 
is usually much larger and may bo as great as 2X10 7 , or 
even more; consequently, flexibility may be thought of in 
terms of parameter A alone. Values of A less than about 
200 usually indicate relatively inflexible rings, and increasing- 
values indicate increasing ring flexibility. 

The right-hand side of equation (1), F(d), contains terms 
resulting from the application of distributed loads; if only 
concentrated loads are present F(0) becomes zero. 


The solution of equation (1) depends upon the nature of 
the roots of the auxiliary algebraic equation associated with 
the differential equation. Of the required eight roots of the 
algebraic equation, four may be real and four complex, or 
all eight may be complex. (For the special case of no 
eccentricity of ring and sheet, sets of roots for several com- 
binations of parameters A and B are given in table I.) 

If there are four real roots ±k u ±k 2 and four complex 
roots ±(a±bi), the solution for the bending moment is 

M= CrfM- f C II e- k ^+C 1 „e^+ 

C y e ta+ii) 0 + C yl e ( °- M) 6 + C v ne ~ la+m) 0 + 

CViiie ~ e J- Particular solution (2) 

which can be written in real form as 

M=C 1 sinh kid-\-C 2 cosh & 1 0+C r 3 sinh k 2 d+Ci cosh k 2 d J r 
(\ sinh ad cos bd-\-(\ sinh ad sin bd+C\ cosh ad sin bd-\- 
C s cosh ad cos bd-\-gJ(6) (3) 

If there arc four pairs of complex roots ±(c 1 ±di'i) and 
±{c 2 ±d 2 i), the bending moment is 

M= C lB +V> 6 + C u e (c i 0 + Cine ~ + "W 9 + C lY e ~ 8 + 
Cye < c 2 -m 2 » »+ C vl e ”+ Cy „e - «+ C vin e ~ ^ "+ 

Particular solution (4) 

which, written in real form, becomes 

M=C X sinh erf cos <7,0-)- C 2 sinh CiO sin di0-\-C 3 cosli c 2 0 sin d\0A 
(\ cosh erf cos di6-\rC s sinh c 2 d cos (U6+ (\ sinh c. 2 d sinrZ 2 0 + 
C 7 cosh c 2 6 sin d 2 d T C s cosh c 2 d cos <7 2 0+>/(0) (5) 

Since the particular solutions gJ(9) and jJ(d) depend upon 
the form of F(0), no formal solutions are given here. Further 


TABLE I.— ROOTS OF AUXILIARY ALGEBRAIC EQUATION WHEN e=0 FOR SEVERAL COMBINATIONS OF A AND B 


\ B 
A 

20 

60 

150 

350 

Root symbols 

10 2 

± 1. 250205 
±7. 745793 

±(0. 612234 ±1. 4679881) 

±0. 956180 
±13. 416404 

±(0. 475554 ±1. 2973971) 

±0. 743842 
±21. 213203 

±(0. 371310 ±1. 1893351) 

±0. 571042 
±32. 403439 

±(0. 290865 ±1. 1131441) 

zb &i 
zb &2 

zb (fl zb bi) 

10 3 

±2. 027436 
± 7. 744229 

±(0. 972538 ±1. 9969011) 

±1. 613689 
±13. 416370 

±(0. 798821 ±1. 7145251) 

±1. 326187 
±21. 213202 

±(0. 661073±1. 5219621) 

±1. 092601 
±32. 403439 

±(0. 547534 ±1. 3740901) 

rfcfci 
zb /C2 

zb (« zb bi) 

10< 

±3. 157270 
±7. 728363 

±(1. 447531 ±2. 8183901) 

±2. 523294 
±13. 416028 

±(1. 236394 ± 2. 3889511) 

±2. 117014 
±21. 213188 

±(1. 052 130 ±2. 0850721) 

±1. 795392 
±32. 403439 

±(0. 896469 ±1. 8457881) 

zb/Ci 
i IC2 

zb (o zb bi) 

10 5 

±5. 016203 
db 7. 539604 

± (2. 058771 ±4. 0301991) 

± 3. 835266 
± 13. 412609 

±(1. 837706 ±3. 4177221) 

± 3. 230542 
±21. 213049 

±(1. 59501 7 ±2. 9595621) 

±2. 769566 
±32. 403431 

±(1. 379274 ±2. 5939991) 

zb 

zb &2 

zb (o rb bi) 


232355—63 2 
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remarks concerning the solution of equation (1) are confined 
to the solution in real form, either equation (3) or equation 
(5). 

Equation (3) or equation (5) expresses the moment in any 
particular panel where A, B, and R are constant; such a 
solution exists for each panel constituting the structure to 
be analyzed. For the case of a circular cylinder, solution 
(3) or (5) is the same as the energy solutions of references 
2 and 3 for a similar two-bay structure. 

RELATION BETWEEN THE MOMENT M AND OTHER FORCES 
AND DISPLACEMENTS 

All forces and displacements in the panel may be expressed 
in terms of the moment given by equation (3) or (5) and 
derivatives of this moment. As previously noted, this 
moment is the moment at the sheet center line, and with the 
exception of the ring bending moment all other forces and 
displacements are those at the center line of the sheet. When 
no eccentricity is involved the sheet center line coincides with 
the neutral axis of the ring. Of major interest in a structure 
such as the one described herein are the forces, stresses 
(obtainable from the forces), and displacements listed in 
the following paragraphs, together with their mathematical 
expressions which are readily obtained in general form in the 
appendix. These expressions become considerably simpli- 
fied in the absence of distributed loads. (For positive sign 
convention of the displacements and forces sec figs. 1 and 
3, respectively.) 



Figure 3— Free-body diagram of ring section showing positive directions of forces. 


The eight quantities listed in this paragraph are associated 
with the boundary conditions (discussed in the next section) . 
These quantities — the forces and displacements at the 
flexible ring, the shear flow at the panel edge, and the axial 
displacement between the flexible and rigid rings at the panel 
edge — are: 

(1) Bending moment, M 

(2) Shear force 


V= 


dM 
'R de 


— m 


(3) Axial force 


T7 1 d 2 M . „ , dm 
• H= ~R -dd^ +Rf+ de 

(4) Sheet shear flow 

1 \dM , d?M T) ,r i ,df , l ( 

q -mUr + ~dF- R V'' + i» + n\. 

(5) Axial displacement 


m- 


d 2 m 

~~d¥ 


)]} 


L 2 (d 2 M , d'M 


r £t\ 


4 R 3 Et'\ de 2 1 dd i 

(6) Tangential displacement 

L {dM d 3 M 
- 2 R 2 Ot\ dd + de 3 

L^ m+ d : M_ R 


T dh d 2 j . 1 / dm d 3 ra\~B 

1 idd + W + Ti\ de + d¥)] { 


A+— + - 
+ de^R 


m- 


dhn 

'ltd 2 


)]} 


w= 


6 R'Et'Xde 3 1 de 5 
(7) Radial displacement 
L 




d 2 h . d 3 f . 1 / d 2 m , d A m 

de 2 + de 3 + R \W + ~d¥ 


+ 

)]} 


{d 2 M , d*M n 2 Vdh . d 2 f , 1 f dm ^d^mW^ ^ 
MBGtXlTe 2 ~ h R\de^ de 3 /Jp 


U \d*M , d % M , J2 r d 3 h , d\f , 1 (d?m , dhn\]\ 
6 R i Et' \ de 3 we 6 11 Idd^ R\d6^ de 5 ) J) 

/ . e d 2 M „ , dm\ 

\ M+ Rle 2 - eRJ - e lie) 


+ 


g(M+ 


(8) Rotational displacements 


L ( dM 


dw v Q _ 
ds + R 2R 3 Gt\ dd 


d 3 M d h M 

W + Tef " 


E, [ k +w+‘M+ 


£j,i 

d e 3 ' R 


, _ d 2 m . d i m\ 

m+2 -W + H¥) 


IJ (d 3 M 


GR 5 Et'\ de 3 


+ 


t d 5 M d‘M 
' de 5 + de r ~ 

, d 4 m d 6 m 

' de 3 + ~d¥ 


R 2 




d 2 h,d i h,d 3 f d\f l / d 2 m 
dJ 2 + de i+ de 3 + de 5 + R \de 2 " 


\"I) e (dM . ed 3 M n df d 2 m 
)\) + El\ do + R de-~- (iR Te- e H¥ 


The shear and normal stresses in the two supporting panels 
follow. The shear stress r in the sheet is given by 


(dM. 

d 3 M 

{ de H 

" de 3 


r 

The stringer normal stress <j is 

-R 2 




L-x(d 2 M . d 4 M 
~ 2 RH'\ de 2 + c W " 


dh . d 2 f . 1 f dm , d 3 mVl) 

Jd + de 2 + Ti\dl + 'dl f )]) U 


The ring bending moment — that is, the bending moment 
about the neutral axis of the ring — is 


, , , , , e d 2 M „ . dm 

M ‘= M +RHiF- eRf - e Te- 


rn 


For no eccentricity the moment at the sheet center line be- 
comes the ring bending moment. 
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BOUNDARY CONDITIONS 

Inasmuch as a general differential equation of the form of 
equation (1) applies for each of the panels composing a 
given structure, a like number of bending-moment solutions 
of the form of equation (3) or (5) results. For the deter- 
mination of the unknown coefficients appearing in these 
solutions a sufficient number of boundary conditions must 
be found, one for each coefficient. Regardless of the number 
of panels, the boundary conditions involve only the eight 
expressions listed in the previous section, consisting of the 
forces and displacements which occur at the sheet center 
line. 

Although, for a particular problem, the application of 
these expressions depends upon the structure to be analyzed 
and the loading, all the boundary conditions may be 
summarized in one general statement: No change in 
displacement or forces can occur across a boundary unless a 
concentrated force is applied at the boundary, in which case 
the change in forces must equal the applied force. Any 
concentrated load in the plane of the ring can be resolved into 
a radial or shearing force, a tangential or normal force, and 
a moment. Then when a concentrated load is applied, the 
boundary conditions require that the difference in shear forces 
of the adjacent panels be equivalent to the applied radial 
load, the difference in axial forces be equivalent to the 
applied tangential load, and the difference in moments equal 
the applied moment. In the absence of any concentrated 
loads, all the forces and displacements must be continuous; 
that is, all eight expressions in one panel must equal the 
corresponding expressions in the adjacent panel. If a con- 
centrated load is applied within a panel, it is necessary in 
the analysis to consider the point of application of the load 
as a boundary and, lienee, to consider the panel as two 
panels, one on either side of the load. 

Although terms associated with distributed loads appear 
in the expressions for the boundary conditions, distributed 
loads affect the boundary conditions only indirectly inas- 
much as they affect the displacements. 

Further discussion of boundary conditions, as they apply 
to one of the numerical examples, is given in the following 
section. 

NUMERICAL EXAMPLES 

The numerical examples deal with two geometrically 
similar cylinders constructed of four panels forming a 
doubly symmetrical, essentially elliptic cross section (fig- 2), 
with the flexible ring subjected to concentrated radial, 
moment, and tangential loads at an intersection with the 
major axis. Cylinder 1 has a very stiff loaded ring and 
cylinder 2, a relatively flexible loaded ring. For each cylinder 
the sheet thickness is constant, there is no eccentricity of 
ring and sheet, and the radius Ri of the top and bottom 
panels is one-third the radius R 2 of the middle panels. The 
moment of inertia of each ring is constant, but because of 
the change in radius, the rings change in relative stiffness 
from one section to another as indicated by the change in 
A/B given in table II. As seen in figure 4, the top panel, 


TABLE II.— NUMERICAL VALUES USED IN COMPUTATIONS 

r 4 _tw p _ Et'Rn 

|_ /L 3 ’ GtL 2 J 


Cylin- 

der 

Section 

A 

B 

A 

B 

Roots of auxiliary 
equation 

1 

1 and 3 

50 

729 

20 

9 

0. 03 

±0.236301 

±2.581786 

± (0. 109747 ±1. 0195321) 

2 

50 

20 

2. 50 

±1.061451 

±7.745880 

±(0.521510± 1.3544931) 


1 and 3 

150000 

729 

20 

9 

92. 59 

± (2.509792 ± 0.621 6651) 
±(0.91 551 3 ±2. 101 9991) 


2 

1 50,000 

20 

7500.00 

±5.534415 
± 7.395633 

±(2.1 81 206 ±4.291 8821) 


e--o° 



Figure 4. — Cross-sectional view of half of loaded ring used for numerical examples. 

section 1, joins the middle panel, section 2, at 0=60°, and 
the middle panel joins the bottom panel, section 3, at 
0=120°. (These dimensions were also used in constructing 
fig- 2.) 

Comparisons of the calculated distributions of bending 
moment, shear force, and axial force in the ring and of 
shear flow in the sheet with the distributions given by the 
elementary theory are shown in figures 5 to 10. The neces- 
sary numerical values used in the calculations are given in 
tables II and III. For ease in reading figures 5 to 10, the 
abscissa, although it actually represents distance along the 
perimeter, is given in degrees measured from the vertical 
axis of symmetry (as shown in fig. 4). Thus, since the ring 
perimeter of section 2 (fig. 4) is three times that of section 1 
or section 3, whereas the angular distance for all three 
sections is the same (60°), the distance in figures 5 to 10 
along the abscissa from 120° to 60° (corresponding to 
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TABLE III.— VALUES OF COEFFICIENTS 1 FOR NUMERICAL EXAMPLES 



Cylinder 1 

Cylinder 2 

Coefficient 

Load 

Load 


Radial 

Moment 

Tangential 

Radial 

Moment 

Tangential 

c. 


0. 147507 

0. 249767 


—0. 860733—03 

0. 119301—03 

c 2 

0. 679031—01 



0. 949137—04 

C 3 


. 275041—04 

— . 157406—03 


. 292731—03 

— . 371518—04 

Ci 

. 125895—05 



— . 284700—03 

C 5 


— . 136156 

. 754809 


— . 148018—01 

. 221238—02 

c 6 

— . 525229 



— . 197121—02 

c 7 


— . 165618 

. 236338—03 


— . 521539—02 

. 484959—03 

c s 

. 282018—01 



— . 354900—02 










C 9 

. 171565—01 

. 188166 

. 202467 

— . 657213—08 

. 140829—06 

— . 161423—07 

C, 0 

. 303490—01 

— . 121785 

— . 111398 

— . 534558—01 

. 270768 

— . 697692—01 

Cu 

— . 264559—10 

. 116591—10 

— . 116591—10 

. 397640—09 

— . 161540—08 

. 710774—09 

C\2 

. 889303 

— . 508263—01 

. 851830 

— . 479259 

— . 380612 + 01 

. 638771 

Cn 

. 329176—01 

— . 296012 

. 472934 

— . 122989—01 

— . 361474—01 

. 485865—02 

Cu 

— . 490646 

. 373871—01 

. 290987 

. 829422—03 

— . 491296—01 

. 702908—02 

c u 

. 192013 

— . 184225 

— . 160403 

— . 117551—02 

. 471184—01 

— . 668234—02 

Cl 6 

. 720834—01 

. 274929 

— . 400188—01 

. 131739—01 

. 392728—01 

— . 527300—02 

C,- 

— . 119460 

— . 164389—01 

. 174960 

- . 119054—05 

. 778086—04 

— . 192398—04 

Cjs 

. 315057 

. 363960 

. 224540 

— . 202194 + 01 

— . 711069 + 01 

. 872951 

Cn 

— . 898224—08 

. 174922—09 

. 301271—07 

. 342801—04 

. 804877—04 

. 168139—04 

^20 

— . 711275—02 

— . 259599 

— . 364881 

. 216297 + 01 

. 882602 + 01 

— . 112279 + 01 

c„ 

— . 640348—01 

. 207820—02 

. 128917+01 

— . 497509—01 

— . 142194 

. 194136—01 

C 22 

— . 531867 

. 131177—01 

. 506058—01 

— . 226239—01 

— . 709217—01 

. 370731—02 

c 23 

— . 254629 

— . 304436 

— . 328941 

. 375148—01 

. 124468 

— . 104012—01 

C M 

. 194401—01 

— . 242335—01 

— . 613555—03 

. 532861—01 

. 167720 

— . 241624—01 


i The following convention is used to indicatelmultiplication factors: +01=X10; +02=X10 2 ; — 01 = X10 -! ; — 02= X10 -2 ; and so forth. 


section 2) represents three times the distance from 180° 
to 120° or from 60° to 0° (corresponding to sections 1 and 3). 

For the numerical examples, the labor of computation 
necessary to calculate the bending moment and other 
desired quantities may be shortened somewhat through cog- 
nizance that antisymmetrieal loading produces an anti- 
symmetrical moment distribution about the vertical axis 
and symmetrical loading yields a symmetrical distribution. 
Hence, only half the cylinder at the ring need bo considered 
(fig. 4). The procedure used in obtaining the numerical 
results is illustrated by taking as an example cylinder 2 
subjected to a concentrated radial load at 0=180°. The 
discussion to follow is confined to cylinder 2 so loaded. 

DIFFERENTIAL EQUATIONS AND SOLUTIONS 

In the top panel, section 1 (fig. 4), the differential equation is 

d*M 14 d 6 M 37 d i M 20 d 2 M 900000 , , 

(/ (9 s 3 de 6 3 de* 3 dd 2 + 729 M_U 


complex roots (see table II) : 

±(ci±dp) = ± (2.509792 ±0.621665i) 

±(c 2 ±rf 2 i) = ± (0.915513 ±2.101999i) 

Since symmetrical loading is applied to the structure, the 
moment in this section is given by only the symmetric terms: 

M iCClion i = C 2 sinh c-fi sin d 1 e+C i cosh ctf cos d]0+ 

C s sinh c 2 0 sin d 2 9+C s cosh c 2 0 cos d 2 d (9) 

In the middle panel, section 2, the differential equation is 


d s M _ 5g «f_ n9 §*_M_ go ^M+9 000 ooM=0 


dd s 


dd* 


dd A 


dd 2 


The auxiliary algebraic equation has the eight roots (see 
table II) : 

±&i = ±5.534415 
±& 2 = ±7.395633 


The resulting auxiliary algebraic equation yields the eight 


±(a±bi)= ±(2.181206±4.291882i) 
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180 120 100 80 60 O 

8, deg 


(a) Ring bending moment. 

(b) Sheet shear flow. 

Figure 9 -Bending-moment and shear-flow distributions produced by concentrated 
tangential load T at 0=180°. 

The moment in section 2 is then 

M, cetlon2 =C 2 sinh hd+C la ' cosh k x d+C n ' sinh k 2 0+ 

C 12 ' cosh k 2 9-\-Cn sinh ad cos 60+ 

Cn sinh ad sin bd+C 15 cosh ad sin 60 + 

C u cosh ad cos bd 

but since k x d and k 2 d, when 0=60° and 0=120°, are large 
enough so that their respective hyperbolic sines and cosines 
are of almost equal magnitude, it is better for computational 
purposes to rewrite the moment in section 2 as 

M, cction 2 =C,e^+C 10 e-^+C n e^+C l2 e-^ + 

Ci 3 sinh ad cos 60+Cu sinh ad sin 60+ 

C\ 5 cosh aO sin 60+CV, cosh ad cos bd (10) 

For the bottom panel, section 3, the differential equation 
is the same as for the top panel and the auxiliary algebraic 
equation has the same roots. Hence, the moment is given by 

M sccticn a= On' sinh c x d cos +0+ C ls ' sinh c x d sin d x d+ 

C n ' cosh c x 9 sin did + CW cosh c x 9 cos d x 9 + 

<7 21 sinh c 2 0 cos d 2 0 + C 22 sinh c 2 d sin d 2 d + 

C 23 cosh c 2 d sin d 2 d +6' 24 cosh c 2 d cos d 2 d 



Figure 10.— Ring shear-force and axial-force distributions produced by concentrated 
tangential load Tat 0 = 180°. 

However, for 0=120° and 0=180° the hyperbolic sines and 
cosines, respectively, of c,0 are almost identical, and it is 
advisable for computational purposes to rewrite the moment 
in the bottom panel in the form 

M accUo ,^—C 17 e c i e cos d x 6 +C 18 e~ c ) B cos di0 + 

(7 19 e c i 6 sin d x 9 +O m e~ c ^ sin d x 9 + 

C n sinh c 2 0 cos d 2 d +C 22 sinh c 2 0 sin d 2 d + 

C 23 cosh c 2 0 sin d 2 d +C 2i cosh c 2 0 cos d 2 d (11) 

BOUNDARY CONDITIONS 

Equations (9) to (11) contain twenty unknown constants; 
hence, twenty boundary conditions arc needed. No bound- 
ary conditions are found at 0=0, since continuity of all 
forces and displacements is already satisfied as a consequence 
of taking advantage of the symmetry of the structure and 
loading. The boundary conditions to be used in the calcu- 
lations must be found at 0=60°, 0=120°, and 0=180°. 

Eight of the boundary conditions are supplied at 0=60°, 
where, in the absence of any concentrated loads, all eight 
expressions in section 1 must equal the corresponding expres- 
sions in section 2. That is, the moment, shear force, dis- 
placements, and so forth must be continuous. For example, 
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the first condition requires continuity of the moments or 

At section 1 Af section 2 0 


which can be written 


n • i C 1 7T • d\ IT • /-> 1 CiTT d\TV . 

C* smh — sin — — b 6 4 cosh — cos —= — b 

do o o 

„ . , Coir ■ doiv , si , Coir do it 

C 6 smh -y- sm — — b C s cosh — - cos — — 


(■ 


k\TT _/Cl7T k?TT _ A'?7T 

C,e s +C 10 e 3 + C n e 3 +C 1 o_e 3 + 

bit 

13 f " 


n • 1 ffilT b IT . ■ 1 d IT ■ V <1 . 

C+ smh — cos -g-+u u smh — sm — + 


„ , (Itt . bir 

Gis cosh — sin — 


d-Cie cosh ^ cos v^)=° 


The solution of the twenty equations given by the boundary 
conditions yields the values of the unknown coefficients C 2 , 
Ci, .. . C 23 , C 24 shown in table III. The bending moment 
in the loaded ring at any angle 0 is then found by using 
equation (9), (10), or (11) with the appropriate coefficients. 
For example, substituting coefficients C, 7 to C 2i from table 
III into equation (11) for 0=180° gives 

ikf 180 .= — 0.119054X10-^^(2. 656444X10 3 )(— 0.372983) + 
(— 0.202194X10 1 PI? 1 )(0.376443X10- 3 )(-0.372983) + 
(0.342801 X lO^PP,) (2.656444 X 10 3 ) (0.927839) + 
(0.21 6297 XIO'PPi) (0.376443 X 10" 3 ) (0.927839) + 
(-0.497509 X 10 _ 1 PPi) (8.444906) (0.949098) + 

( — 0.226239 XI O-’PPO (8.444906) (0.3 14983) + 
(0.375148 XlO-'PPj) (8.901257) (0.314983) + 

(0.532861 X lO-'PP,) (8.901257) (0.949098) 

= 0.161390PP 1 


The other seven boundary conditions appearing at 0=60° 
are written in a similar manner. The same continuity of 
the forces and displacements must exist at 0=120°. Thus 
sixteen boundary conditions are provided, with the four 
remaining conditions to be found at 0=180°. 

For an applied radial load, the sum of the shearing forces 
on either side of the cut must equal the applied shearing 
load P. Because of the symmetry of the structure and load- 
ing, half of this load is carried by section 3 and the other 
half of the load is carried by the panel to the left, which 
need not be considered. Furthermore, the shear flow in the 
sheet is 0 at 0=180°, there can be no tangential or sidewise 
displacement at 0=180°, and there are no angular displace- 
ments caused by a concentrated radial load at this point. 
Hence, the four boundary conditions needed may be summed 


up as follows : 

V i8o» — 0.5 P 

(12a) 


O' 180” — 0 

(12b) 



(12c) 


(t+7rL=° 

(12 cl) 


Equations (12a) 
substitution to 


to ( 12 d) can be reduced by proper 



=0.5Pi?i 
= -0.5 PR 
=0.5 PIP 
= —0.5 PR 


i 


For example, equation ( 12 a) can be written 

C i 7 e c ' \c 1 cos d i 7 r — d 1 sin d ity — C^e ~ c 1 T (ci cos d x ir + d ! sin d 1 + + 

C'i 9 e ,:i ’ r (c 1 sin d^ + di cos diir ) — C 20 e~ c ' ir (ci sin di-w—d, cos d] 7 r) + 

C 2 \{c 2 cosh c 2 ir cos doir — d 2 sinh c 2 ir sin rf 2 7 c)~b 

C 22 (c 2 cosh c 2 ir sin d 2 ir+<f 2 sinh c 2 -ir cos d 2 ir) + 

Cwfe sinh c 2 t sin do_irAd 2 cosh c 2 ir cos d 2 P) + 

C 2 i (c 2 sinh c 2 7 r cos d 2 ir — d 2 cosh c 2 x sin rf 2 7 r) = 0.5PP 1 


CONCLUDING REMARKS 


The results obtained from the numerical examples agree 
with those previously obtained for circular cylinders (see, 
for example, references 2 and 3) in indicating that concen- 
trated loads applied to flexible rings produce stresses in the 
rings and shell considerably different from those computed 
from an engineering analysis (wherein the ring is treated as 
a free ring supported by the usual elementary torsion and 
bending shears). Ring flexibility is essentially indicated by 


/ y/;«\ 

the parameter + ( = -j-ff ) 


; where A is everywhere less than 


about 200 the engineering analysis is adequate, but if A 
exceeds 200 such an analysis is inadequate. 

The main effects of flexibility are to change considerably 
the distribution of stress and the magnitudes of the maximum 
stresses. The change in maximum stresses is indicated in 
the following table, which gives the approximate ratios of 
the absolute values of the maximum stress coefficients for 
cylinder 2 (A — 206 near the load but 150,000 some distance 
away) to the maximum stress coefficients obtained from an 
engineering analysis : 


N \ Type of 
\ stress 

Load 

Cm 

C, 

tv 

On 

Radial 

0.5 

4.5 

1.0 

1.4 

Moment 

1.0 

47.5 

3.4 

fi. 3 

Tangential 

.2 

5.4 


1.0 


For much larger values of A — that is, for greater flexibility — 
the ratios greater than unity would increase considerably, 
those less than unity would decrease somewhat, and the 
ratios of unity would remain unchanged. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., August 1, 1951. 




APPENDIX 

DEVELOPMENT OF GENERAL DIFFERENTIAL EQUATION 


Consider the two-bay panel with uniform geometrical 
properties (fig. 1), loaded in the plane of the middle ring. 
In accordance with the assumption that the shear stress in 
the sheet may vary circumferentially but remains constant 
in the axial direction, an infinitesimal element of the sheet 
is subject to the stresses shown in figure 1. 

Equilibrium in the x-direction requires that 


^-tdsdx-\-^t' dxds = 0 
os ox 


(Al) 


Since 

and 

equation (Al) becomes 


rt=q 

at'—p 


*2 + *l = o 

dx^ds 


(A2) 


Integrating equation (A2) with respect to x gives 

However, since the end ring supports are free to warp out 
of their planes, p = 0 at x=L, and therefore 


A(-)=i §f 


Then 


V=(L — r) 


A M 


ds 


The strain relation in the ^-direction gives 
d u 

so that 


dx 


JL 

Et' 


d u 1_ __ . d q 

d x~Et' [ X) ds 

Integrating equation (A3) with respect to x gives 


(A3) 


-w{ L 


x 2 \ d q 
2 / ds 


+M*) 


but since, from symmetry of the structure and loading, 
u= 0 at x=0, 

/ 2 (s) = 0 

so that 


u= 




2 7 ds 


The shear strain can be expressed as 

_ q _d« . dr 
^ Gt ds dx 


(A4) 


(A5) 


from which u can be eliminated and v can be found in terms 
of q. Differentiating equation (A4) with respect to s and 


substituting the result into equation (A5) gives 

L ( Lr x ‘ 2 '\ b2( L 

d x~Gt Et'\ X 2 ) ds 2 

Integrating with respect to x yields 

q 1 (Lx 2 x 3 \ b 2 q , , , , 
v ~Gt x Et’ \ 2 QJ ds 2+ ^ 3< - 

but since the end ring supports are rigid in their planes, 
r=0 at x=L, and 


, , , L , L° d 2 (z 
J&S)— 2 


3 Et' ds 2 


so that 


v — ~Gi q QEt'^ Lx2 ** 2i3) ds ?2 


This displacement in the sheet at the ring, where x=0, 
becomes 

L , 1 L 3 d 2 q 


v 0 - Gt q_- 


\\R 2 EV M 2 


(AO) 


The relation between the strain in the sheet in the circum- 
ferential direction and the deflections at the ring is 


_l (dv _ o 

\de 


-w 


where by continuity w is the radial deflection of both ring 
and sheet. Since the strain e s in the sheet is the strain at the 
outermost fibers of the ring, 

o’., e M r 

fs ~E~~ El 

Thus the radial deflection can be expressed as 

eRM R . dv 0 

w= -ET- + de 


(A7) 


Differentiating equation (AO) and substituting the result into 
equation (A7) yields 


w = 


eRM a L dq, 1? d 3 q 


(AS) 


El Gt dd 1 liREt' do 3 

The relation between the ring bending moment and the 
radial deflection is given by (see, for example, reference 11) 


El ( d 2 w . \ 


(A9) 


Performing the indicated differentiation of equation (A8) 
and substituting the result into equation (A9) eliminates w 
and gives 

El VeR (d 2 M R 


f El VeR (d 2 M R ,, \ L / d 3 q dq\ + 
Mr (R-ef [_EI V do 1 +jX1r ) Gt\dd 3+ dd) + 


1 ? (d 5 q . d 3 q\ 
3R 2 Et' \dd 3 'd,6 3 )_ 


(A10) 
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It is now convenient to find the relation between the shear 
flow in the sheet and the moment at the sheet center line, 
as well as to express the moment at the ring center line in 
terms of the moment at the sheet center line. Figure 3 
shows an infinitesimal segment of the loaded ring and the 
positive directions of the forces. For the first of the three 
equations of equilibrium, taking moments about the origin 
gives 

— dM-\-RdH-\-R 2 h dd+Rrn dd- J r2R 2 q dd— 0 
or 

(A11) 

Summation of the tangential forces (neglecting terms of 
higher order) gives 

V dO—Rh dd~2Rq dd—dll— 0 
or 

JTJ 

v = / r j \-Rh \ 2Rq (A 12) 

Summation of the forces in the radial direction gives 

— H dd-{-Rj dd — dV=0 
or 

H=-^+Rj (A13) 


Substituting the value of from equation (Al2) into equa- 
tion (All) yields the shear force at d 


ldM 
R de 


m 


(A 14) 


Differentiating equation (Al4) and substituting the result 
into equation (Al3) gives the value of the normal force at 
the sheet center line at 8 


H= 


1 d 2 M 
R de 2 


Rf+ 


dm 

W 


(A 15) 


Differentiating equation (A15) and substituting the result 
into equation (All) gives the expression for the shear flow 
in the sheet at 9 


1 \dM , d 3 M r>2 Ti, i df , 1 / i d 2 m\~\) , A ,„ S 

2 2 R 2 ( de + de 3 R L h+ de + R ( m+ de 2 (A16) 

The bending moment about the ring center line can be seen 
to be composed of the moment M and the product of the 
axial force H and the eccentricity e or 


Mr— M—eH 


, f , e d 2 M „ . dm 


(A17) 


Differentiation of equations (Al6) and (A17) and sub- 
stitution of the result into equation (AlO) leads to the 
desired differential equation 


d s M 
d9 8 


(2-35) 


d s M 
de 6 ' 


■(> -« b +t) 


d A M 
de 4 


o r„ 2eA( e e 2 \]d 2 M . . ( e e 2 \ 

3 L j? r \ 2 ~r + ~r 2 )]i^ + 6A \ 1 ~r + w) m ~ F( '^ 


where 


(A 18) 


r W =B‘[g+a-3B)g-3sg]+ 


R 2 


R 




d e f 
de 6 

6 eA 

R 

( d 7 m 


+ (1-3 B) 


dM 

de 4 


o / d 2e 2 A\ d\f 
3 \ B 5 2 7d+ + 


(i — D + 


e e 
R^R 2 




+ 


6e 2 +\ d 3 m 


|^ +(2 _ 3B ,^ + ( 1 - 6£+ ^) 


de 3 


~R + R 2 )\~d6) 


Equation (Al8) is the general form of the differential equa- 
tion for the bending moment at the sheet center line for any 
of the panels composing the structure. The bending 
moment in the ring for any such panel is given by equation 
(A17). 
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